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ENERGY DEPENDENCE OF HORIZONTALLY 
INCIDENT ATMOSPHERIC MUONS 


Kaichi Maeda 


ABSTRACT 

Horizontally incident atmospheric muons and their parent particles travel 
one of the longest paths in the earth's atmosphere, leading to the largest possi- 
ble atmospheric effects on the muon production process. Diffusion equations for 
cosmic ray mesons in the atmosphere are inadequate to treat these effects at 
large zenithal direction. The calculation is, therefore, made by modifying the 
Chapman function for the non-isothermal atmosphere, and computations are also 
applied to a mountain altitude, specifically for the laboratory at Mt. Chacaltaya 
(5200 m). 
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ENERGY DEPENDENCE OF HORIZONTALLY 
INCIDENT ATMOSPHERIC MUONS 


1. INTRODUCTION 

The purposes of this paper are: first, to extend the author's previous calcu- 
lation (Maeda 19G4)* to higher energy and larger zenith angles; and second, to 
apply these extended calculations to the determination of the effects of atmos- 
pheric structure on the zenith angle distribution of the atm ^spheric muons. In 
particular, the direction of che maximum muon intensity as a function of muor. 
energies at two different altitudes will be discussed. The atmospheric muon is 
defined by Wolfendale (1968). 

2. COSMIC RAYS IN THE ATMOSPHERE 

In M-l, variations of cosmic ray intensity in the atmosphere are expressed 
in integral form as a function of the atmospheric depth, x in g/cm 2 instead of 
the solutions of the so-called diffusion equations of cosmic rays in the atmosphere. 
In this expression, the following three features— which the diffusion equation can- 
not describe properly— are taken into account rigorously: (i) variation of air 
density with altitude, corresponding to changes in the vertical thermal structure 
of the atmosphere, (ii) variation of zenith angle along the straight oblique path 
with respect to the local horizontal due to the curvature of the earth's surface 
(see Figure 1), and (iii) a non-linear energy-range relation in the high energy 
regions due to the radiation loss and the direct pair-production loss, which 
dominate at extremely high energies. As was shown in M-l, these integrals 
agree with several approximate solutions of the so-called diffusion equations. 

In other words, the solution of these different equations corresponds to some 
special cases of the present results under several approximations. Notations 
and definitions in this paper are mostly the same as those used in M-l. 


2. 1 Differential Intensity of Cosmic Ray Mesons 

The intensity of cosmic ray mesons, n. ( E it x, 6 *) with energy between E ( 
and E i + dE ( at the depth x with zenith angle between 6* and 6* + dO* is given 


“This will be referred to as M-l, later in this report. 


by 
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sec 9 * ( x ' ) dx ' 


( 2 . 1 ) 


where 




sec 6 * ( x " ) 


dx" 
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E. w (x")p(x"X 


sec 0* ( x" ) dx" 


The suffix i stands for pion ( tt ) and kaon (k), respectively, and F^V) is the 
differential production spectrum of cosmic -ray mesons. The energy E /' at the 
depth x" (at the production level) is related to Ej at the depth x by a range- 
energy relation. 

2.2 Energy-Relation Between Muons and Par ent Mesons 

As was shown in M-l, the differential energy soectrum n (E (i , x, 9 *) of 
muons at production is given by 


f r i B sec 9 * 

>V(V Xl<? *) = 'e P(E„)n i (E, , x. 9 ") dE, (2.2) 

\ 

where 



(2.3) 


2.3 Differential Intensity of the Atmospheric Muons 

The differential energy spectrum i (E, x Q , 9 ) of muons at the depth x Q , 
arriving from the direction with zenith angle 9 , is different from that of the 
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production n (E , x, 9 *) at the depth x, but is connected by the so-called " sur- 
vival probability," W(E, x Q , x,<9), in the following way: 




x, 9 *) W(E, x Ql x, 6 ) dx , 


( 2 . 0 ) 


where the local zenith angle 9 * (x) at the depth x along the straight path of the 
muon, arriving at the depth x Q with the zenith angle is related to , as can 
be seen from Figure 1, by 


sin 0 * (x) 
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(2.7) 


The survival probability W(E, x 0 , x, 6 ) is given by 


W (E, x q , x, 9 ) 
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0 sec 9 * ( x ' ) 
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( 2 . 8 ) 


2.4 The Energy Loss in the Atmosphere 

Detailed studies of energy loss processes of muons have been made by sev- 
eral cosmic ray workers up to around 1000 GeV (Barrett et al., 1952; Ashton, 
1961; Ozaki, 1962; Dayman et al., 1963). These works are, however, based on 
the underground experiments and results are expressed with respect to the 
Standard Rock. Since the energy loss process of muons has been known as es- 
sentially the same as that of electrons but with larger mass, their results, which 
depend on Z/A and Z 2 /A of the medium, can be reduced to those in air by using 
Table 1. In Table 1 the density p, atomic number Z, atomic weight A, Z/A and 
Z 2 /A of air and of the Standard Rock are shown for comparison. 

The change in a muon’s energy along its penetration path in the atmosphere 
can be calculated by using the average range-energy relation in air, which is 
derived from the empirical formula for the energy loss rate of muons, i.e., 


where df dx see 0* (x) is a path element of the muon in the atmosphere, and 
the constants a e , b p are 2,5 x 10~ 3 (GeV/g cm 2 ) and 2,5 * 10“ 6 (g/cm 2 )' 1 , re- 
spectively (M-l). 


3. HORIZONTALLY INCIDENT COSMIC RAYS 

At the mountain laboratories, the cosmic rays arriving from below the 
horizontal direction can be observed. 

To calculate the intensities of these cosmic rays, the expression described 
in the previous section is not suitable because of the sec 6* term, which diverges 
at 6 = 90° as the altitude h(x) approaches that of the observation point, z. The 
calculation is made, therefore, by using the Chapman function instead of sec 6* . 
The Chapman function is defined by the constant scale height atmosphere, i.e., 
the isothermal atmosphere (Chapman, 1931). This is, however, easily extended 
for the standard atmosphere, in which the isothermal stratosphere is situated 
above the polytropic (constant lapse rate) troposphere. The atmospheric depth, 
x (in g cm’ 2 ) in the standard atmosphere is shown in Figure 2 as a function of 
altitude, z (in km). 

3.1 In the Isothermal Atmosphere 

All expressions in the previous sections are expressed by the Chapman 
function, the integral of ch(0, k) with respect to the local zenith angle, k t along 
a straight path of muons in the atmosphere as follows: 
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and 


E U (M ar E + (a e + b e E)p 0 ch(0 t k^dk 1 . 


(3.7) 


3.2 In the Standard Atmosphere 

Because of the difference in the atmospheric depth-height relation between 
the upper and lower atmosphere, i.e., above and below h = h t , all terms con- 
taining ch(i9, \) in the expressions in Section 3.1 must be replaced by ch(b, k) 
for\<\ t , andc^(a, k) for respectively. Corresponding to Equation 

(3.2), therefore, we get 


n^E.W) = F(E/) -Vj ' e ' D s (E ‘ A ' ch(b, V) 
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(3.12) 


Similarly, the relation between the muon’s energy E (X) at the production level 
h( x ) and E at the observing altitude z with zenithal angle, 9 is, instead of Equa- 
tion (3.7), given by 


E u (\) = E + (a e + b e E) 


•X re 


'X 


ch(b, X ' ) dX' + A. 


c %( a , X ' ) dX ' 


• (3.13) 


It should be noted that ch(b, X) andc^(a, X) are defined only for X <X t and 
Q — ^ — respectively, i.e., 

ch(b, X) = 0 for X > X t , 

c4(a, X) = 0 for X < X t . (3.14) 
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It should also be noted that most data from muon intensity measurements are 
integral spectra which are given by 




x Q , d) • dE 


(3.15) 


where E 0 is the cut-off energy defined by conditions of the muon-detector and 
V (E, x o> ) is given by (3.1). 

4. NUMERICAL CALCULATIONS AND RESULTS 

4.1 In the Isothermal Atmosphere 

Directional differential intensities of the atmospheric muons, i (E, x Q , ), 
which are calculated by using Equations (3.1) -- (3.7), are shown in Figure 3(a) 
for x Q = 1030 g/cm 2 (z - 0 km, sea level) and (b) for x Q = 530 g/cm 2 (z = 5.2 
km). In these figures, intensities are plotted (in GeV~ 1 cm’ 2 sec” 1 ster” *) 
against the zenith angle, , with energy E (in GeV) at the observing point as a 
parameter. In all cases, the attenuation mean free path L n of the primary cos- 
mic rays and the nuclear absorption mean free path of L n pions are both as- 
sumed to be 120 g/cm 2 . To see the difference between warm and cold atmos- 
pheric conditions, results are also shown for two temperatures, T = 340° and 
T = 170°K as heavy and light curves, respectively. The warm and cold atmos- 
pheres correspond to scale heights of 10 km and 5 km, respectively. 

As can be seen from these figures, the horizontally incident atmospheric 
muon intensity in the isothermal atmosphere decreases with the atmospheric 
temperature, both at sea level and at mountain altitudes. It should be noted, 
however, that in near-vertical directions, the high energy muon intensity in- 
creases with atmospheric temperature increase; this lias been known as the 
’’positive temperature effect" on the cosmic ray intensity (Duperier, 1951; 

Barrett et al., 1952). 

On the other hand, at low energies (below around 100 GeV), the muon intensity 
in near-vertical directions decreases with increase of atmospheric temperature, 
and this is well known as the "negative temperature effect" on cosmic ray intensity 
(Blackett, 1938; Maeda and Wada, 1954; Dorman, 1957). 
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In Figure 4, the energy dependence of directional differential muon intensity, 
* (E, x 0 , 0) is shown (a) for sea level and (b) for z = 5.2 km. Corresponding 

to Figures 3(a) and (b), two cases of the atmospheric temperature, 340°K and 
170°K arc shown by light and heavy lines, respectively, in Figures 4(a) and (b). 
Since the energy dependence of i (E, x 0 , ) is indicated in Figure 3, only two 

cases of the zenith angle are shown in Figure 4, i.e„, (a) for sea level 0 = 0°, 

85°, and 90°; and (b) for z = 5.2 km altitude, = 0°, 90°, and 92.3°, respectively. 

One conspicuous feature of i fL (E, x Q , 8) is the direction of the maximum 
muon intensity, which is not only a function of energy, E, but also of atmospheric 
depth (i.e., altitude) x Q , of temperature T, of the atmospheric structure, and of 
the mean free path L i of the muons, parent particles. This is shown in Figure 5, 
where the zenith angle of the maximum muon intensity, max (E), is plotted 
against the muon energy at the observation point, E (in GeV), for sea level and 
mountain altitudes (5.2 km). To see the dependence on the atmospheric temper- 
ature, two cases, T = 340°K (H = 10 km) and T = 170°K (H = 5 km), are pre- 
sented. For a comparison, one case for the standard atmosphere is also shown 
in Figure 5. From this figure, one can see that: (i) The direction of the max- 
imum muon intensity is almost independent of energy at very high energies 
(above around 10 TeV). (ii) At low energies (below around 100 GeV), the direc- 
tion of the maximum approaches zenith rapidly with decreasing energy, (iii) The 
zenith angle of the maximum muon intensity is larger at mountain altitudes than 
at sea level for the same muon energy at the observing point, (iv) The zenith 
angle of maximum intensity is also larger in the cold atmosphere than in the 
warm atmosphere, (v) Obviously, the dependence of $ max on the atmospheric 
structure is more sensitive in the lower energies than in the higher; and in the 
latter, # m ix approaches 90° almost independently of energy E. 

4.2 In the Standard Atmosphere 


Because of the existence of a vertical temperature gradient in the lower part 
of the actual earth’s atmosphere, the rate of increase of atmospheric depth, x, with 
decreasing altitude, h, is less in the lower atmosphere (the troposphere) than in 
the upper atmosphere (the isothermal stratosphere), as is shown in Figure 2. 

Since the decay rate of muons in flight is inversely proportional to their 
energy, which decreases with penetration depth in the atmosphere, the intensity 
of muons arriving deep in the atmosphere with a given initial energy differs 
from that calculated for the isothermal atmosphere, especially in the near- 
horizontal directions. 

Corresponding to Figure 3, the directional integral intensity of atmospheric 
muons produced by pions I (E 0 , x Q , 6) is shown in Figure 6(a) for sea level and 
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(b) for a mountain altitude, ^ = 5.2 km. Similarly, in Figure 7, I (E 0 , x Q , r ) 
corresponding to kaon-produced muons it: shown (a) for sea level and (b) for a 
mountain altitude, z = 5.2 km. 

Figure 8 shows the energy dependence of the zenith angle of the maximum 
muon intensity, ' max (E) at (a) sea level and (b) a mountain altitude z = 5.2 km. 
To see the difference in # max due to different values of the nuclear absorption 
mean free path of muon-producing mesons, three eases are indicated, assuming 
L t to be equal to 1/2 L n , L n and 2L n where Lj is the attenuation mean free path 
of pion-producing cosmic-ray primaries, i.e., i = -tt for pions and i - k for 
kaons. Figure 8 indicates that the smaller L j} the larger (E) for the same 
energy E. In other words, the zenith angle dependence of atmospheric muons of 
a given energy produced by hypothetical mesons with a small absorption mean 
free path (i.e., a larger cross section for the nuclear absorption) is the same 
as that of the muons produced by the same mesons with a larger mean free path 
(i.e., small cross section for nuclear absorption in air) at higher energy. Since 
these zenith angle dependences of the atmospheric muons are due to the com- 
petition process between the nuclear absorption of muon-producing mesons, 
which is proportional to l/L ( , and the decay in flight of these mesons in the up- 
per atmosphere, which is proportional to EF/Ex, the result for small L f can be 
reproduced approximately by large L i with larger energy, as is discussed in 
M-l. The decay factor of the kaon is larger than that of the pion. Thus, to 
produce the same zenith angle dependence of muon intensity, the energy of the 
kaon must be roughly seven times larger than that of the pion. This is also 
shown in Figure 8 assuming L k = 1/2 L n , L k - L n , and L k = 2L n . 


5. DISCUSSION 

An important matter related to the present calculation is the so-called 
"sec 0-law” for the zenithal distribution of cosmic-ray muons (Bergeson et al., 
1967; 1968). As can be seen from Figures 3-6 this is only approximately cor- 
rect for small zenith angles. 

The difference in the zenith angle distribution of high energy muons between 
sea level and mountain altitudes (z = 5.2 km) shown in Figures 3-6, can be sum- 
marized as follows: (i) Below around 500 GeV, the zenithal distribution, par- 
ticularly the direction of the maximum muon intensity, changes with the absorp- 
tion mean free path of the muons’ parent particles as well as with the abundance 
of kaons vs pions. (ii) Above 1 TeV, the zenith angle dependence and the direc- 
tion of the maximum muon intensity are almost independent of the characteristics 
of parent particles and rather dependent on the atmospheric structure. 

Further details of the present calculations on the atmospheric muons will be 
discussed elsewhere (Maeda, 1970). 
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Table 1 


Physical Constants of the Air and the Standard Rock 


Average Value 

The Air 

The Standard Rock 

p, density 

1.25 x 10' 3 (g/cm 3 ) 

2.65 (g/cm 3 ) 

Z, atomic number 

7.37** 

11. 

A, atomic weight 

14.78** 

22. 

Z/A* 

0.499 

0.5 

Z 2 /A* 

3.675 

5.5 


* The average values of Z/A andZ /A are calculated by X. fjZj /A £ and respectively, 

where f j is fractional weight. 

"Rossi (1952), p. 295. 
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Figure 2. Comparison of the atmospheric depth-altitude relation in the standard atmosphere and 
empirical formula. Dashed lines correspond to the extrapolations of the formula beyond appli- 
cable altitude range. 
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Figure 3. Zenith angle dependence of the differential intensity of atmospheric muons in the iso- 
thermal atmosphere, i^(E, x Q , 6 ) in cm 2 sec 1 ster 1 GeV \ plotted against zenith angle 6 , 
wi th the energy of muons at the observation point, E (in GeV), as a parameter, (a) for sea level, 
Xg = 1030 g/cm 2 (z - 0 km) and (b) for mountain altitudes Xg - 530 g/cm (z - 5.2 km). Heavy 
and light curves correspond to the intensity in the warm (T = 340 K, H “ 10 km)and cold atmos- 
phere (T = 170°K, H - 5 km), respectively. 
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Figure 6. Zenith angle dependence of the directional intensity of atmospheric muons 
(inteqral spectra) produced by the decay of pions in the Standard Atmosphere 1^ (Eq, 
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g/cm^). The cut-off energy of the spectrum is E q. 
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Figure 8(a). Zenith angle of the maximum atmospheric muon intensity in the standard atmosphere, max , p.otted 
against to the energy of muons at the observation point for sea •evel (z 0). L^and L k are taken as a parameter 
of each curve, where heavy and light curves represent the pion-produced and kaon-produced muons, respectively. 
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Figure 8(b). Zenith angle of the maximum atmospheric muon intensity in the standard atmosphere, max , plotted 
against to the energy of muons at the observation point for mountain altitudes (z 5.2 km). L^and are taken 
as a parameter of each curve, where heavy and light curves represent the pi;«n-prod uced and kaon-produced 
muons, respectively. 




